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1.  INTRODUCTION 


This  research  program  has  been  directed  towards  obtaining  a fundamental 
understanding  of  chemical  kinetic 'gas  dynamic/part icle  interactions  in  rocket 
nozzle  and  exhaust  plume  flows.  The  work  is  motivated  primarily  by  the  need 
for  accurate  predictions  of  such  rocket  plume  observables  as  IR  radiation, 
radar  cross-section  and  visible  smoke.  The  nozzle  flow  studies  have  emphasized 
the  determination  of  gas  and  particle  properties  at  the  exit  plane  which  are 
used  as  input  to  rocket  plume  codes;  the  plume  work  has  primarily  concentrated 
on  determining  IR  radiation  properties  of  low  altitude  afterburning  plumes  from 
subscale  liquid  propellant  motors.  An  additional  effort  has  been  directed 
toward  the  calculation  of  particle  flows  in  afterburning  plumes  with  the  ulti- 
mate goal  of  predicting  the  visibility  of  low  signature  propellant  plumes. 


II.  ACCOMPLISHMENTS 


A.  Fully-Coupled  Rocket  Nozzle  Flow  Model /Computer  Code 

A significant  advance  in  nozzle  flow  calculations  was  achieved  with 
the  development,  under  this  contract,  of  a fully-coupled  rocket  nozzle  flow 
model  and  computer  code  (FULLNOZ) . 1 The  model  simultaneously  treats  chemical 
kinetics  and  partlcle/gas  thermal  and  dynamic  nonequilibrium  effects  for  con- 
toured nozzles,  including  the  influence  of  the  wall  boundary  layer.  The  com- 
puter code  is  based  on  the  MULTITUBE  code  developed  by  Boynton,2  which  utilizes 
the  streamtube  method  (as  described  by  Boynton  and  Thomson*)  to  integrate  the 
hyperbolic  governing  equations  of  steady  supersonic  flow.  Basically,  this 
method  solves  the  corresponding  finite-difference  equations  along  and  perpendic- 
ular to  streamlines.  The  elliptic  Navier-Stokes  equations  are  reduced  to 
hyperbolic  form  by  assuming  that  diffusional  effects  along  streamlines  are 
small  compared  to  diffusion  across  streamlines.  This  is  a very  good  assumption 
for  nozzle  flows  and  enables  one  to  solve  an  initial  value  problem  (where  a 
marching  procedure  can  be  used),  rather  than  the  more  difficult  boundary  value 
problem. 


The  operation  of  FULLNOZ  is  relatively  straightforward,  starting  with 
the  specification  of  initial  gas  and  particle  properties  just  downstream  of  the 
throat  (Mach  number  > 1)  and  the  nozzle  wall  contour.  Additional  input  data 


1.  Pergaraent,  H.S.  and  Thorpe,  R.D.,  "A  Computer  Code  for  Fully-Coupled 
Rocket  Nozzle  Flows  (FULLNOZ),”  AeroChem  TP-322,  AFOSR  TR-75-1563, 

April  1975. 

2.  Boynton,  F.P.,  "The  MULTITUBE  Supersonic  Flow  Computer  Code,"  General 
Dynamics/Convair  GDC-DBB  67-003,  February  1967. 

3.  Boynton,  F.P.  and  Thomson,  A.,  "Numerical  Computation,  of  Steady,  Super- 
sonic, Two-Dimensional  Gas  Flow  in  Natural  Coordinates,"  J.  Computational 
Phys.  3,  379-398  (1969). 


2 


include  the  thermodynamic  properties  of  all  constituent  cheaical  species  (which 
are  Incorporated  via  curve  fits  of  JAKXAF  data),  a choaical  reaction  awe  ha  n Isa 
and  rate  coefficients,  the  physical  properties  of  the  particles  and  the  nozzle 
wall  temperature  (for  the  boundary  layer  calculations'.  The  marching  scheme 
proceeds  downztreaa,  froa  one  orthogonal  (to  the  flow  atreamllne)  surface  to 
the  next,  computing  flow  properties  and  coaposltion  within  a specified  nuaber 
of  streaatubes.  At  each  orthogonal  surface  the  code  also  calculates  wall  shoar 
stress  and  heat  transfer,  boundary  layer  displacement  thickness  and  velocity 
and  temperature  profiles,  coupling  the  shear  stress  and  heat  transfer  to  the 
'wall'  at reamt ube  properties.  The  turbulent  boundary  layer  analysis  utilizes 
the  Van  Driest  transformations4  and  the  experlaental  data  of  Keener  and  Hopkins,* 
which  relate  the  compressible  skin  friction  coefficient  to  aeasured  velocity/ 
temperature  profiles  in  flows  with  favorable  pressure  gradients.  Coupling  be- 
tween the  boundary  layer  and  'wall'  streamtube  properties  is  achieved  via  solu- 
tion of  the  momentum  Integral  equation  to  give  boundary  layer  momentum  thickness 
as  a function  of  distance  along  the  nozzle  wall.  The  displacement  thickness  is 
then  determined  from  the  velocity  profile  and  the  momentum  thickness. 

FULLNOZ  has  now  been  established  as  an  operational  code  for  analyzing 
rocket  nozzle  flows.  It  has  been  delivered  to  the  Air  Force  Rocket  Propulsion 
Laboratory  (AFRPL)  and,  via  AFRPL.  is  available  throughout  the  country  to  Indus- 
trial and  government  laboratories.  It  will  be  used  by  Lockheed/Palo  Alto  as 
part  of  their  visible  plume  signature  study  for  AFRPL.  AeroC-hem  has  used  it 
extensively  for  AFRPL  (under  subcontract  to  Aerodyne  Research,  Inc.)  to  predict 
nozzle  exit  plane  properties  for  a number  of  solid  and  liquid  propellant  tacti- 
cal and  strategic  missiles.  These  properties  were  then  used  by  AeroChem  as  in- 
put to  the  calculation  of  exhaust  plume  electrical  properties  (electron  density 
and  electron-neutral  collision  frequency)  which  in  turn  were  input  to  plume 
radar  cross  section  calculations  made  by  Aerodyne.  In  addition  we  have  used 
FULLNOZ  to  calculate  gas  and  particle  properties  within  the  Space  Shuttle  solid 
rocket  booster  nozzle.* 

The  speed  and  versatility  of  FULLNOZ  are  demonstrated  in  Table  I, 
which  summarizes  input  conditions  and  computer  run  times  for  several  types  of 
systems  of  interest  to  the  Air  Force.  The  usefulness  of  FULLNOZ  can  be  further 


4.  Van  Driest,  F..R.  , "Turbulent  Boundary  Layer  in  Compressible  Fluids," 

J.  Aeron.  Sci.  18,  1A5-160  (1951). 

5.  Keener,  E.R.  and  Hopkins,  E.J.,  "Van  Driest  Generalization  Applied  to 
Turbulent  Skin  Friction  and  Velocity  Profiles  Measured  on  the  Wall  of  a 
Mach  7.4  Wind  Tunnel,"  AIAA  J.  U,  1784-1785  (1973). 

6.  Pergament,  H.S.,  Thorpe,  R.D.,  and  Hwang,  B.,  "N0X  Deposited  in  the 
Stratosphere  by  the  Space  Shuctle  Solid  Rocket  Motors,"  Final  Summary 
Report,  Phase  II,  AeroChem  TP-333,  NASA  CR-144928,  December  1975;  see  also, 
Thorpe,  R.D.,  Pergament,  H.S.  and  Hwang,  B.,  "N0X  Deposition  in  the 
Stratosphere  by  the  S ;ace  Shuttle  Solid  Rocket  Motors,"  JANNAF  9th  Plume 
Technology  Meeting,  CPIA  Publ.  277  (Applied  Physics  Lab.,  Johns  Hopkins 
Univ.,  Silver  Spring,  April  1976),  pp.  317-352. 
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enhanced  by  Ch«  lncorpotatlon  of  a transonic  flov  solution  <e.g. . r'let  of  Nicker- 
•on  and  Klle<>elr  to  allow  the  code  to  start  at  the  combustion  cheaper  rst'ier 
than  downstream  of  the  throat.  Other  modifications  that  can  easily  be  made  v’ich 
would  extend  the  capability  of  FULLNOZ  Include:  (1'  low  Reynolds  number 

flows,  where  the  boundary  layer  is  a large  portion  of  the  nozzle  flow,  (2'  noz- 
zle film  cooling,  l.e.,  the  effects  of  mass  injection  on  wall  heat  transfer  and 
shear  stress,  and  (3)  nozzle-induced  shocks. 

B.  Interpretation  of  IR  Radiation  Data  on  Subscsle  Liquid  Propellant  Motor 

Plumes 


Data  obtained  by  the  Air  Force  Armament  Laboratory*  on  subscale 
kerosene /oxygen  rocket  motors  (fired  into  a vacuum  tank)  have  shown  that  mea- 
sured IR  plume  radiation  In  the  4-5. 5u  band  is  in  sharp  disagreement  with  pre- 
dictions using  state-of-the-art  modeling  techniques.  These  techniques  Include: 
(1)  thermochemical  equilibrium  chamber  and  nozzle  calculations,  (2)  a coaxial 
turbulent  mlr.lng/aonequilibrlua  chemistry  exhaust  plume  analysis  and  (")  IR 
radiation  ’band  model’  calculations  with  CO,  as  the  radlstlng  species.  The 


7.  Nickerson,  G.R.  and  Kllegel,  J.R.,  "Axisymmetrlc  Two-Phase  Perfect  Gas 
Performance  Program,"  TRW  Systems  Report  No.  02374-6006-R000,  Vols.  I and 
II,  April  1967. 

8.  Ebeoglu,  D.B.  and  Martin,  C.W.,  "Experimental  Verification  of  Infrared 
Plume  Predictions  for  a Rocket  Engine,"  AFATL-TR-74-191,  November  1974; 
see  also,  Ebeoglu,  D.B.,  Martin,  C.U.,  Brail,  S.,  Victor,  A.,  Mantz,  J., 
and  Rothschild,  W.J.,  AIAA  Paper  No.  75-1231;  AIAA/SAE  11th  Propulsion 
Conference,  Anaheim,  California,  29  Sept.  - 1 Oct.  1975. 
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nature  of  t he  discrepancy  between  measured  end  predicted  radiation  Is  that  for 
0/F  ratios  less  thau  stoichiometric  the  total  Tdlant  Intensity  is  such  greater 
than  that  predicted  by  tho  model,  while  for  high  o/F  ratios  the  model  over- 
orsdicts  the  data.  The  goal  of  our  work  under  this  contract  war  to  explain  the 
discrepancies  between  theory  and  experiment  and  to  develop  a nodal  which  can 
reliably  be  used  to  pr«>dict  the  plune  gas  dynamic  and  chemical  properties. 

A careful  evaluation  of  all  assumptions  utilixed  in  the  original  pre- 
dictions demonstrated  that  the  tost  probable  cause  for  the  failure  of  the  model 
is  that  it  did  not  account  for  the  effects  of  the  (observed)  Incomplete  combus- 
tion in  the  rocket  motor  chamber.  Based  on  this  observation  a new  'phenomeno- 
logical' model  was  devsloped  which  assumes  that  some  of  the  fuel  does  not  burn 
in  the  chamber  but  la  pvrolyxed  to  a lower  order  hydrocarbon  before  leaving  the 
rocket  noxxle.  The  amount  of  fual  unturned  is  then  determined  by  comparing  the 
measured  rocket  performance  (the  C*  efficiency)  to  thet  predicted  via  thermo- 
chemical  equilibrium  calculations.  This  determination  is  made  by  assuming: 

(1)  s fixed  percentage  of  the  fuel  does  not  burn  but  forms  a specific  pyrolysis 
product  (e.g.,  C.Hi.,  C,Hfc  or  CH*,  etc.'  from  the  liu  lid  kerosene  and  (2s  the 
pyrolysis  product  remains  inert  in  the  chamber  and  noxxle  expansion.  The  per- 
centage of  urburned  fuel  leaving  the  noxxle  is  assumed  to  be  that  amount  for 
which  the  measured  and  calculated  C*  efficiencies  are  equal.  The  model  further 
assumes  that,  upon  entering  the  plume,  the  unburned  fuel  participates  in  three 
overall  kinetic  processes:  gas  phase  oxidation,  polymerisation  to  form  soot 

and  soot  oxidation.  These  global  processes  are  coupled  to  the  conventional  ten- 
step  CO/Ha  oxidation  scheme  through  the  following  three  reactions. 

Gas-phase  Oxidation: 


(CHa)  + OH  - CO  + H3  + H 


Soot  Formation: 


n(CHa)  - ( CH  g ) 

n 


Soot  Oxidations 


(CHa)n  + n02  * nCO  + nHaO 

where  (C!ia)  represents  an  arbitrary  fuel  fragment. 

The  results  obtained  to  date*  ^shovn  in  Fig.  1)  have  incorporated 
only  the  gas-phase  oxidation  step.  For  0/F  ■ 1.8  ‘ (O/F^stoich  * 3*5)  the 
results  Indicate  that  account!  .g  for  the  effects  of  unburned  fuel  vi>»  the  ?N?ve 
model  will  give  substantially  increased  plume  tampemturas  and  IR  radiation  if 
the  turbulent  mixing  rate  (ambient  ai**  with  exl.-»us<-  nroducts)  is  sue’.'  that  all 


9.  P^rgament,  H.S.,  Thorpe,  R.D.,  and  Miller,  W.J.,  "interpretation  of  IR 
Radiation  Plume  Data  on  Subscale  Kerosene/Oa  Rocket  Motors,"  JANNAF  9th 
Plume  Technology  Meeting,  CPIA  Publ.  277  (Applied  Physics  Lab..  Johns 
Hopkins  Univ.,  Silver  Spring,  April  1976),  pp.  125-143. 
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Fig.  1 Comparison  between  measured  station  mediation  (O'T  - 1.") 
and  prediction?  rith  phenomenological  model. 

Original  predictions  from  Ref.  3;  c » eddy  viscosity 

coefficient;  (CHa)  + OH  **  CO  + H?  + H; 

* 1 x 10“11  exp(-7000/RT) ; (CHa)  is  unburned  fuel. 

the  'unburned  fuel'  burns  in  the  plume.  This  increased  level  of  IR  radiation 
is  in  much  better  agreement  with  the  IR  data  than  the  original  predictions, 
demonstrating  the  general  validity  of  the  new  phenomenological  model. 

Much  work  remains  to  be  done  before  the  new  model  can  be  generally 
accepted  as  a viable  calculational  procedure  for  this  type  of  afterburning  plume. 
This  work  should  include: 

(1)  Comparisons  between  predictions  and  data  over  the  entire  range  of  0/F 
for  which  radiacion  measurements  were  made. 

(2)  Establishing  an  eddy  viscosity  coeffici  mt  consistent  with  the  actual 
plume  velocity  field. 

(3)  Determining  the  influence  of  inaccuracies  in  measured  Cgff  on  the  pre- 
dictions, via  a series  of  parametric  calculations  in  which  the  % un 
burned  fuel  is  varied  over  a reasonable  range. 
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(4)  Incorporating  the  aoot  formation  and  oxidation  steps  in  tha  sfterburninj: 
pluat  raac t ion  mechanism. 

(5)  Compar iaooa  between  tlaorr  and  tiparlatar.  for  tha  aora  recect  aat  of 

data  obtained  at  APATL.***  Tbaaa  data  ware  taker  on  tha  same  not  or  a, 
but  with  naaa  flows  of  5 g/sec  (in  contrast  to  9 g/sec  for  tha  tast 
data  Interpreted  in  the  present  work)  and  burned  barttane  and 

■ethane  (CHfc),  in  addition  to  kerosene. 

Because  un burned  hydrocarbons  entering  the  pliae  are  apparently 
primarily  responsible  for  the  large  Ik  signatures  of  these  pluses  for  rich  mix- 
tures. the  following  practical  cooclualon  can  be  drawn:  The  'rocket  motor* 

used  for  the  simulation  studies  should  be  designed  a priori  aa  a hydrocarbon 
fuel  injector,  with  (possibly)  direct  injection  into  the  plume  rather  thin  with 
unburned  fuel  injected  as  a result  of  Incomplete  combustion  in  the  chr  : «r . 

C . Particls  Motion  in  Afterburning  Rocket  Plumes 

The  purpose  of  this  task  was  to  develop  a computational  tool  to 
describe  particle  motion  in  afterburning  rocket  plumes.  Specifically,  we  have 
Incorporated  the  particle  subroutine  in  FULLNOZ  (subroutine  DRAG),  which  numer- 
ically evaluates  the  gas/part icle  interaction  terms.  Into  the  AeroChem  Low 
Altitude  Plume  Program'*  (LAPP),  the  standard  tool  throughout  the  country  for 
computing  afterburning  plume  properties.  The  gas  phasa  momentum  and  enargy 
equations  in  LAPP  have  Seen  modified  to  account  for  these  Interaction  terms  and 
the  particle  conservation  equations  have  been  formulated  in  streamline  coordi- 
nates for  use  in  Lapp.  These  equations  must  be  coded  and  the  LAPP  input  format 
modified  to  accept  particle  data  in  order  to  develop  an  operational  code.  Once 
these  tasks  have  been  completed,  the  LAPF  code  will  be  able  to  accept  input 
data  on  nozzle  exit  plane  gas  and  particle  properties  and  routinely  compute 
gas  and  particle  temperatures  and  velocities  in  afterburning  plumes. 


T These  comparisons  are  being  made  as  part  of  the  Cooperative  Plume  Modeling 
Program  under  the  auspices  of  the  TTCP  W-4  panel.  Further  information  may 
be  obtained  from  Capt.  W.J.  Rothschild,  DLMQ,  Air  Force  Armament  Laboratory, 
Eglin  AFB,  FL. 

10.  Rothschild,  W.J.  and  Martin,  C.W.,  "Experimevital  Results  for  Hydrocarbon 
Exhaust  Infrared  Model  Verification,"  AFATL-TR-76-74,  July  1976. 

11.  Mikatarian,  R.R. , Kau,  C.J.,  and  “ergament,  H.S.,  "A  Fast  Computer  Program 
for  Nonequilibrium  Rocket  Plume  Predictions,"  Final  Report,  AeroChem  TP-282, 
AFRPL-TR-7 2-94 , NT1S  AD  751  984,  August  1972. 
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